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H 2 emission in high-z ultra-luminous infrared galaxies 
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ABSTRACT 

Using the IRAM Plateau de Bure interferometer (PdBI), we report the detection of water vapor in six new lensed ultra- 
luminous starburst galaxies at high redshift, discovered in the Herschel Astrophysical Terahertz Large Area Survey 
(H- ATLAS). The sources are detected either in the 2o2-ln or 2n-2o2 H2O emission lines with integrated line fluxes 
ranging from 1.8 to 14 Jy km s _1 . The corresponding apparent luminosities are /iLn 2 o ~ 3 — 12xlO 8 L0, where \i is 
the lensing magnification factor (3 < fi < 12). These results confirm that H2O lines are among the strongest molecular 
lines in high-z ultra-luminous starburst galaxies, with intensities almost comparable to those of the high-J CO lines, 
and similar profiles and line widths (^200-900 km s _1 ). With the current sensitivity of the PdBI, the water lines can 
therefore easily be detected in high-z lensed galaxies (with F(500/im) > 100 mJy) discovered in the Herschel surveys. 
Correcting the luminosities for amplification, using existing lensing models, Lh 2 o is found to have a strong dependence 
on the infrared luminosity, varying as ^Lj^ 2 . This relation, which needs to be confirmed with better statistics, may 
indicate a role of radiative (infrared) excitation of the H2O lines, and implies that high- z galaxies with Lir > 10 13 Lq 
tend to be very strong emitters in water vapor, that have no equivalent in the local universe. 
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Water plays an important role in the warm dense inter- 
stellar medium of galaxies. First, after CO, H2O is the 
most abundant oxygen-bearing molecule, and, second, it 
can be an important coolant of the warm gas. Due to the 
Earth's atmosphere, bulk gas-phase water can only be de- 
tected from space or from the ground toward distant objects 
with redshifts that move the H2O lines into atmospheric 
windows. 

The study of water emission lines in nearby galaxies has 
recently made significant progress thanks to the availability 
of their infrared/submillimeter spectra using the spectrom- 
eter mode of the Herschel Spectral and Photometric Imag- 
ing REceiver (SPIRE, Griffin et al. 2010) and the Photode- 
tecting Array Camera and Spectrometer (PACS, Poglitsch 
et al. 2010). The spectra of local ultra-luminous infrared 
galaxies (ULIRGs) and composite AGN/starbursts such as 
Mrk231 (van der Werf et al. 2010), Arp220 (Rangwala 
et al. 2011; Gonzalez- Alfonso et al. 2012), and NGC4418 
(Gonzalez- Alfonso et al. 2012) reveal a wealth of water 
lines and the presence of associated molecules such as OH + , 



* Herschel (Pilbratt et al. (2010) is an ESA space observatory 
with science instruments provided by European-led Principal 
Investigator consortia and with important participation from 
NASA. 



H20 + and isotopologues. Together with high-J CO lines, 
these lines provide an important diagnostic of the warm 
dense cores of nearby ULIRGs. 

In the cases of Mrk231 and Arp220, water emission 
lines up to energy levels of E up /k = 642 K are detected with 
strong line fluxes that reach 25%-75% of the neighboring 
CO emission line fluxes. Spectral surveys made with Her- 
schel show that \ow-z ULIRGs always exhibit bright H2O 
lines, whereas, only one third of the sample of luminous 
infrared galaxies (LIRGs) displays luminous H2O emission 
lines (van der Werf et al. in preparation [prep.]), indicat- 
ing that the strength of the water lines and the infrared 
luminosity, Lir, must be related. The analysis of the H2O 
emission lines in Mrk231 shows that the excitation of the 
water molecules results from a combination of collisions and 
infrared excitation through far-infrared lines in warm dense 
gas (> 100 K, > 10 5 cm -3 ). Moreover, the far-infrared ra- 
diation field dominates the excitation of the high levels and 
their emission lines (Gonzalez- Alfonso et al. 2010). Pre- 
liminary evidence from the comparison of the spectra of 
Arp220, Mrk231 and NGC4418 shows that the properties 
of the water emission lines in their nuclear regions vary as 
a function of chemistry, nucleosynthesis, and inner motions 
(outflow/infall) - see Rangwala et al. (2011) and Gonzalez- 
Alfonso et al. (2012). 
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First detections of H2O megamasers at high redshift 
were reported in objects at z = 0.66 and 2.64 by Bar- 
vainis & Antonucci (2005) and Impellizzeri et al. (2008), 
while further searches have failed since then (McKean et 
al. 2011). Following initial attempts to detect H2O rota- 
tional emission lines from high-z galaxies (Wagg et al. 2006; 
Riechers et al. 2006, 2009) and tentative detections in IRAS 
F10214+4724 at z = 2.23 (Casoli et al. 1994) and in the 
Cloverleaf at z = 2.56 (Bradford et al. 2009), a series of 
robust detections of non-maser H2O emission lines were re- 
cently reported in high-z sources. The sources detected in 
H2O emission lines include a strongly lensed galaxy, HAT- 
LAS J090302.9-014127 (SDP.17b) at z = 2.3 (Omont et 
al. 2011), which was uncovered in the Herschel Astrophys- 
ical Terahertz Large Area Survey (H- ATLAS) (Eales et al. 
2010); the gravitationally lensed quasar APM 08279+5255 
at z = 3.9 (van der Werf et al. 2011; Bradford et al. 
2011; Lis et al. 2011); HLSJ091828.6+514223, a lensed 
submillimeter galaxy at z = 5.2 in the field of Abell 773 
(Combes et al. 2012), which was found in the Herschel 
Lensing Survey (Egami et al. 2010), IRAS F10214+4724 
(Riechers et al., in prep.). For most of these sources, only 
one water emission line was reported, with the exception 
of APM 08279+5255, where a total of at least five emission 
lines with E u /k = 101 to 454 K were detected (van der Werf 
et al. 2011; Lis et al. 2011; Bradford et al. 2011). These re- 
sults underline the unique and powerful diagnostic power 
of H2O emission lines, which give better insight into lo- 
cal conditions in distant galaxies than may be obtained by 
other means. As in Mrk 231, they reveal the presence of ex- 
tended, warm and dense gas located in the infrared-opaque 
regions of the galactic cores. 

This new window in the exploration of high-z sources is 
based on the combined availability of new instrumentation 
with improved sensitivities at key facilities and the increas- 
ing number of gravitationally lensed sources discovered in 
the Herschel and South Pole Telescope (SPT) cosmological 
surveys (see, e.g., Negrello et al. 2010; Vieira et al. 2010). 
Following the results reported by Omont et al. (2011), we 
present here a new study of water emission in six high-z 
lensed ULIRGs, which were selected from the H- ATLAS 
survey. The data, which were obtained using the IRAM 
Plateau de Bure Interferometer (PdBI), clearly show water 
emission lines in all the sources. Based on these results, we 
derive a clear relation between the H2O and the infrared 
luminosities in high-z ULIRGs. 

Throughout this paper, we adopt a cosmology with 
H = 71kms- 1 Mpc- 1 , Q M = 0.27, n A = 0.73 (Spergel 
et al. 2003). 

2. Sample selection and observations 

The detection of strong H2O emission in the lensed H- 
ATLAS galaxy SDP.17b by Omont et al. (2011) suggested 
that the H-ATLAS survey offered a unique opportunity 
to select an homogeneous sample of bright lensed galax- 
ies spanning a broad range in luminosity and redshift to 
further study the properties of water emission in the early 
universe. We selected from the H-ATLAS catalog, sources 
with strong flux densities (with Fsoo/im > 150 mJy) that 
were well characterized, i.e. had CO redshift measurements, 
additional sub/millimeter imaging and available deflector 
identification (see, e.g., Negrello et al. 2010; Bussmann et 
al. in prep.). In the selection, we also somewhat preferred 



sources at z > 3 to shift the H2O lines in the low frequency 
bands that are easier to observe. Finally, we chose sources 
spanning a wide range in intrinsic infrared luminosities from 

- 5 x 10 12 to - 2 x 1O 13 L . 

In this paper, we report the first results of this survey 
and describe the properties of the H2O emission for six 
new sources. The sample includes (see Table 1 for IAU 
names of the sources) i) sources reported in the initial H- 
ATLAS Science Demonstration Phase paper by Negrello et 
al. (2010) - SDP.9 and SDP.81, together with SDP.17b, 
which was discussed in Omont et al. (2011); ii) two other 
well-studied sources from the H-ATLAS equatorial fields 

- G12.v2.30 (Fu et al. 2012) and G15.v2.779 (Cox et al. 
2011; Bussmann et al. 2012); hi) two sources from the H- 
ATLAS NGP field, NA.vl.144 and NB.vl.78, for which CO 
observations (Harris et al. 2012; Riechers et al. in prep.; and 
this paper) and submillimeter imaging (Bussmann et al. in 
prep.) are available. Table 1 provides details on the H2O 
observations of the seven H-ATLAS lensed galaxies, Table 2 
lists their submillimeter and infrared properties together 
with the estimated amplification factors (see Sects. 3.2 and 
4.1), and Table 3 gives their CO properties available from 
recent measurements using the PdBI. 

The redshifts of the seven selected galaxies range from 
1.58 to 4.24, with the majority in the range 2.0 < z < 3.3. 
In the sample, G15.v2.779 is the only source at z > 4 
and, for the time being, the highest redshift lensed galaxy 
spectroscoptically confirmed in the H-ATLAS survey. For 
each source, at least one of the two strongest low-excitation 
H2O lines is in an atmospheric window observable with 
the PdBI, either H 2 O(2ii-2 02 ) with ^ rest =752.033 GHz and 
E up = 137K, or H 2 O(2 02 -lii) with ^ rest =987.927 GHz and 
Eup = 101 K. As both lines have comparable intensities in 
Arp220 and Mrk 231 (see Sec. 4), in the rare cases where 
both lines are observable, we chose the line at the lower 
frequency that is easier to observe, except in the case of 
SDP.17b (Omont et al. 2011), where the line H 2 O(2 02 -lii) 
was selected to confirm the tentative detection that was re- 
ported by Lupu et al. (2012). Note that these two observed 
lines are both para-H20 and they are adjacent lines in the 
H2O level diagram. 

The H2O observations were conducted in the compact 
D-configuration from December 2011 to March 2012 in con- 
ditions of good atmospheric phase stability (seeing of 1") 
and reasonable transparency (PWV < 1mm). Except for 
SDP.81, which was observed with five antennas, all other 
sources were observed with the six antennas of the PdBI 
array. The total on-source integration ranges from ~1 to 2 
hours (Table 1). We used Bands 2, 3 and 4, which cover 
the frequency ranges 129-174 GHz, 201-267 GHz and 277- 
371 GHz, and the band centers were tuned to the redshifted 
frequency of the selected H2O emission line using the red- 
shifts estimated from CO observations (Table 3). The cor- 
relator (WideX) provided a contiguous frequency coverage 
of 3.6 GHz in dual polarization with a fixed channel spacing 
of 1.95 MHz, allowing us to detect the continuum as well 
as any additional emission lines, if present. In the com- 
pact configuration, the baselines extend from 24 to 179 m, 
resulting in synthesized beams of ~ 1.5" x 1.0" to ~ 4" x 
3" (Table 1). Only SDP.81, G12.v2.30 and NB.vl.78 are 
resolved at this angular resolution (Fig. 2). 

During the observations, the phase and bandpass were 
calibrated by measuring standards sources that are regu- 
larly monitored at the IRAM PdBI, including 3C84, 3C279, 
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Fig. 1. Spectra of the H2O and CO emission lines observed at PdBI toward the seven lensed high-redshift ultra-luminous 
galaxies discussed in this paper. The H2O emission lines are displayed in the left column (showing the (2n — 202) or the (2o2 — In) 
transition; the transitions are identified in each box. The CO spectra are displayed in the right column (showing 3-2, 4-3, 5-4 and 
7-6 rotational transitions - for G15.v2.779, the [CI] emission line is included in the band). The transitions are indentifled in each 
box. All spectra include the continuum, and the H2O frequencies corresponding to the zero velocity are given in Table 4. The red 
lines trace the Gaussian fits to the H2O and CO line profiles and the fit to the continuum emission level. In all cases, the match 
between the H2O and CO spectra is excellent both in redshift and in the details of the line profile - see text. 



MWC349, CRL618, and 0851+202. The accuracy of the 
flux calibration is estimated to range from ~10% at 2 mm 
to -20% at 0.8 mm. 

To complement the H2O data, we include in this pa- 
per new results on the CO emission for the seven sources. 
For five of them, the CO data were obtained in 2011 and 
2012 using the PdBI in the A-configuration. Those data 
were acquired in the frame of a survey of lensed ULIRGs to 
map their CO and dust continuum emission at high-angular 
resolution (Cox, Ivison et al. in prep.) and a full descrip- 
tion of the results will be given in that paper. For two of 
the sources, NB.vl.78 and G12.v2.30, the CO data were 
obtained in August 2012 in the D-configuration with five 
antennas. The observations were made in good weather 
conditions and more details are provided in Table 3. In 
this paper, we only present the global CO spectra (i.e., in- 
tegrated over the source's extent) with the goal to compare 
the characteristics of the H2O and CO emission line profiles. 
A detailed discussion of the morphology, the dynamics, and 
the lensing of these sources will be given in Cox, Ivison et 
al. (in prep.). 

3. Properties of the H2O emission lines 

3.1. General properties 

The seven lensed high-z ULIRGs are all detected with high 
signal-to-noise ratios (S/N > 6, except SDP.81, either in the 
2o2 _ ln or the 2n-2o2 H2O emission line (see Fig. 1), as well 
as in the underlying redshifted submillimeter continuum 
emission (S/N > 20). The water emission lines are strong, 
with integrated fluxes ranging from 1.8 to 14 Jy km s _1 
(Fig. 1 and Table 4). This indicates that in high-z ULIRGs, 
they are the strongest molecular submillimeter emission 
lines after those of high- J CO (see Sec. 4.3; note, however, 
that 2 mm and 3 mm CO lines displayed in Fig. 1 may have 
weaker intensity than the 1 mm H2O lines). 

The H 2 linewidths (FWHM) range from 140 to 
900kms _1 (Table 4). The lines have a variety of profiles 
including single Gaussian profiles (NA.vl.144, NB.vl.78), 
double-peaked profiles (G12.v2.779), or asymmetrical pro- 
files (SDP.9, SDP.17b, SDP.81). The H 2 and CO line 
profiles share similar properties (shape, linewidth) in all 
sources, suggesting that there is no strong differential lens- 
ing effects between the CO and the H2O emission lines. 

From the spatially integrated H2O line flux, Ih20? the 
apparent H2O luminosity, /iLH 2 o where \i is the lensing 
magnification factor, can be derived using the relations 
given in, e.g., Solomon, Downes & Radford (1992). /iLH 2 o 
varies by only a factor of ~ 4 from SDP.81 to NB.vl.78 or 
G12.v2.30 (Table 4). Fig. 3 displays the relation between 
/iLH 2 o and //Lir, including the values for local ULIRGs 
(Yang et al. in prep.). The infrared and H2O apparent 
luminosities are one or two orders of magnitude higher in 
these high-z galaxies than in local ULIRGs. However, the 
ratios Lh 2 o /Lir are similar, although slightly higher for the 
high-z galaxies. 

Figure 2 displays the images of the H2O emission line 
and the corresponding submillimeter continuum emission 
for the seven lensed high-z ULIRGs. The relatively low 
angular resolution (at best ~ 1") of the current data lim- 
its for most of the sources a detailed study of the spa- 
tial properties of their signal. Four of the sources re- 
main unresolved at the present angular resolution: SDP.9, 



SDP.17b, NA.vl.144, and G15.v2.779. The three remaining 
sources are extended and display distinct lensed morpholo- 
gies including two-image configuration system (SDP.81), a 
complex elongated structure (G12.v2.30), and an extended 
structure (NB.vl.78). Each source is described in detail 
below. 

Table 4 reports the values of the H2O velocity- integrated 
flux density (line flux) detected at the peak of the source in 
one synthesized beam, lH 2 opk, and the velocity-integrated 
line flux integrated over the source size, Ih 2 o- In most 
cases, the ratio lH 2 o/lH 2 opk is lower than 1.5, except for 
the extended sources G12.v2.30 and NB.vl.78. For the 
other sources studied in this paper, except for SDP.81, it 
is thus unlikely that significant flux is missed. This is also 
confirmed when comparing our total spatially integrated 
continuum flux density S u with measurements published 
elsewhere at the nearby frequency of 250 GHz, using the 
bolometer MAMBO camera at the IRAM 30-meter tele- 
scope (Negrello et al. 2010; Dannerbauer et al. in prep.). 
NA.vl.144 and SDP.81, which were observed at the PdBI 
close to 250 GHz, agree very well with MAMBO observa- 
tions. For SDP.9 and SDP.17b, which were observed at the 
PdBI close to 300 GHz, the comparison is more difficult; 
however, is somewhat higher than the values expected 
from extrapolations of the MAMBO values. Although there 
is no MAMBO observation for NB.vl.78, agrees fairly 
well with the flux density measured at 880 /im using the 
Smithsonian Millimeter Array (SMA) (Bussmann et al. in 
prep.). Finally, the continuum flux density measured at 
143 GHz for G15.v2.779 is in excellent agreement with the 
value reported by Cox et al. (2011) at 154 GHz. 

The 3.6 GHz bandwidth covered by the PdBI receivers 
(and correlator), which was generally centered on the red- 
shifted frequency of the H2O emission line, encompasses 
other potential emission lines of interest, including i) 13 CO, 
H^ 8 0, H^ 7 0, H 3 0+ in the vicinity of H 2 O(2 02 -lii); and 
ii) H 2 O + (2 2-lii,J3/2,3/2 and J 5 /2,3/2) in the vicinity of 
H20(2n-2o2))- However, the current sensitivity of our 
spectra is too low to detect these lines, which are not seen 
in the Herschel spectra of Mrk231 and Arp220 except for 
the two H20 + lines, which are present in the spectrum of 
Arp220 (Rangwala et al. 2011). Apart from the H2O emis- 
sion lines, no other emission line is detected at the S/N 
ratio of the spectra of the seven high-z galaxies studied in 
this paper. 

3.2. Individual sources 

The following describes the properties of the individual 
lensed galaxies. 

3.2.1. SDP.81 at z = 3.040 

The source SDP.81 is at the highest redshift among the 
SDP lensed galaxies of Negrello et al. (2010). It displays 
two arcs in the SMA 880 /im data, which are also seen in the 
HST (Negrello et al. in prep.) and the PdBI CO imaging 
results (Cox, Ivison et al. in prep.). Preliminary lensing 
modeling indicated a very high lensing magnification factor 
of /i = 19 (Negrello et al. 2010), but this value has been 
revised downward to /1 = 9.5=b0.5 based on recent iKTdata 
(Dye et al. in prep.). This is confirmed by the SMA 880 /im 
imaging results which yield /i = 9.5 ± 0.8 (Bussmann et 
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Fig. 2. Images of the submillimeter continuum (left panel) and H2O line emission (right panel - see Fig. 1 and Table 4 for the 
identification of the transitions) for the seven high-z lensed ultra- luminous galaxies studied in this paper (sources are identified in 
the left panels). Synthesized beams are shown in the lower left corner of each map. The map centers correspond to the coordinates 
given in Table 1 and are indicated by a cross. For the continuum maps, the contours start from 3a in steps of 3a for SDP.81 and 
G12.v2.30, and in steps of 6a for the other sources; for the H2O line emission maps, the contours start from 3a in steps of la (for 
SDP.81) and 2a for the other sources. For each source, the la continuum (mJy/beam) and emission line (Jy km s _1 /beam) noise 
levels are as follows: SDP.81 (0.60 / 0.18), NA.vl.144 (0.35 / 0.31), SDP.9 (0.50 / 0.59), NB.vl.78 (0.28 / 0.30), G12.v2.30 (0.28 
/ 0.28), G15.v2.779 (0.11 / 0.24), and SDP.17b (0.73 / 0.56). Four sources remain unresolved at the current angular resolution, 
whereas SDP.81, NB.vl.78, and G12.v2.30 display complex lensing morphologies both in the continuum and line emission that are 
revealed even at the present low angular resolution. Detailed comments on the sources are given in the text (Section 3). 



al. in prep. - hereafter B13). As discussed in Sec. 4.1, 
we adopted the values of \i derived from the 880 /im SMA 
data by B13 as basic reference for the lensing magnification 
factors listed in Table 2. 

The CO line profile is asymmetrical with a red com- 
ponent that is much stronger than the blue one. The line 



profile of the H2O emission line is also dominated by a main 
red component (peaking at ~ +200 km s -1 with a width of 
Av = 140 km s -1 ). Due to the low S/N ratio of the H 2 
spectrum, the blue emission component of the profile (peak- 
ing at — 100 km s -1 ) that is clearly seen in the CO spectrum 
is only barely visible in the H2O emission spectrum (Fig. 1). 



Table 1. Observation Log 



IAU Name 


ID 


RA a 


Dec a 


Date 


frequency 


Beam(") 


ton(h) C 


HATLAS J090311. 6+003906 


SDP.81 


09:03:11.61 


+00:39:06.7 


2011 


Dec 


244.890 


2.8x1.8 


1.2 


HATLAS J133649. 9+291801 


NA.vl.144 


13:36:50.00 


+29:17:59.6 


2012 


Mar 


235.330 


1.6x1.1 


0.8 


HATLAS J090740. 0-004200 


SDP.9 


09:07:40.05 


-00:41:59.5 


2011 


Dec 


292.550 


3.3x1.6 


1.4 


HATLAS J133008. 4+245900 


NB.vl.78 


13:30:08.56 


+24:58:58.3 


2012 Jan & Mar 


240.802 


1.4x1.0 


2.3 


HATLAS J114637.9-001132 


G12.v2.30 


11:46:37.99 


-00:11:32.0 


2011 


Mar 


232.734 


2.0x1.0 


1.2 


HATLAS J142413. 9+022303 


G15.v2.779 


14:24:13.90 


+02:23:04.0 


2011 


Dec 


143.940 


4.2x3.1 


1.7 


HATLAS J090302. 9-014127 


SDP.17b* 


09:03:03.02 


-01:41:26.9 


2011 


Jan 


298.148 


2.9x1.1 


0.6 



a: J2000 coordinates of the centers of the maps displayed in Fig. 2. 
b: Central observed frequency (GHz). 

c: Total on-source integration time for the PdBI array with six antennas. 
*: The observations of SDP.17b are reported in Omont et al. (2011). 



Table 2. Submillimeter and infrared properties of the lensed ultra-luminous galaxies 



Source ID 


F250 


F350 


F500 






Lir 




(mJy) 


(mJy) 


(mJy) 


(1O 13 L ) 




(1O 12 L ) 


SDP.81 


129±20 


182±28 


166±27 


5.8 


9.5 ±0.8 


6.1 


NA.vl.144 


295±45 


294±45 


191±31 


5.7 


5.3±2.9 


11 


SDP.9 


485±73 


323±49 


175±28 


4.4 


8.5 ± 1.8 


5.2 


NB.vl.78 


273±42 


282±43 


214±33 


11 


10.5 ± 1.4 


10 


G12.v2.30 


323±49 


378±57 


298±45 


15.7 


9.6 ±0.9 


16 


G15.v2.779 


115±19 


192±30 


204±32 


8.5 


4.1 ±0.2 


21 


SDP.17b 


328±50 


308±47 


220±34 


6.9 


4.3±1.2 


16 



Note: F250, F350, and F500 are the SPIRE flux densities at 250, 350, and 500 /im; /iLiR is the apparent total infrared luminosity 
(8-1000 /im); /i is the adopted lensing magnification factor derived from lensing modeling of 880 /im continuum, mostly from 
Bussmann et al. in prep, (see text Sec. 3 for details); Lir is the intrinsic infrared luminosity, derived from the literature (Harris et 
al. 2012; Lupu et al. 2012; Fu et al. 2012; Bussmann et al. 2012) except for NB.vl.78, for which an approximate value is inferred 
from scaling the value of G12.v2.30 with the ratio of the 350 /im flux densities. 

Table 3. Observed parameters of the CO emission lines 



Source ID 


z 


CO line 


Ico 


AVco 


Other CO data 








(Jy.km/s) 


(km s -1 ) 




SDP.81 


3.040 


5-4 


7.0±0.4 


560±40 


Harris et al. (2012) 


NA.vl.144 


2.202 


4-3 


12.3±0.6 


210±10 


Harris et al. (2012) 


SDP.9 


1.574 


3-2 


9.2±0.3 


420±10 


Iono et al. (2012) 


NB.vl.78 


3.111 


3-2 


3.3±0.5 


560±70 


Riechers et al. (in prep.) 


G12.v2.30 


3.259 


3-2 


7.8±1.7 


650±160 


Fu et al. (2012) 


G15.v2.779 


4.244 


7-6 


7.5± 0.6 


840±50 


Cox et al. (2011) 


SDP.17b 


2.305 


4-3 


9.1±0.3 


320±10 


Harris et al. (2012) 



Note: z is the redshift measured from the PdBI CO spectra (Fig. 1); Ico is the integrated flux of the CO line; AV co is the CO 
linewidth (FWHM); references to other CO measurements from the literature are given in the last column. 



Despite the low angular resolution, the dust continuum 
image is clearly resolved (Fig. 2), showing both gravita- 
tional arcs of the 880 /im image of Negrello et al. (2010), 
with a stronger emission in the eastern arc. This struc- 
ture is also seen in the CO map of Cox, Ivison et al. which 
indicates that the stronger red emission traces the bright 
eastern arc, whereas the weaker blue emission corresponds 



to the western arc. Due to the low S/N ratio of the data, 
the H2O emission is seen only in the eastern arc, consis- 
tent with the only clear detection of the red part of the 
spectrum. 

The H2O line flux Ih20 and the apparent Lh 2 o luminos- 
ity, which are given in Table 4, take into account the total 
emission including the red and blue components. However, 



Table 4. Observed parameters of the H2O emission lines and continuum emission 



Source 


Line a 


b 

^obs 


S,pk c 




S^H 2 opk d 


AV H 20 


I H2 opk d 


T d 
1 H 2 


/iL H2 o e 






[GHz] 


[mJy/beam] 


[mJy] 


[mJy/beam] 


[km s _1 ] 


[Jy km s _1 /beam] 


[Jy km s- 1 ] 


[1O 7 L ] 


SDP.81 


2 


244.5000 


12.5 


27.0±1.2 


10.7±1.9 


140±50 


1.6 


1.8±0.5 


32±9 


NA.vl.144 


1 


234.8343 


9.4 


16.4±0.7 


22.8±1.7 


200±50 


4.9 


7.5±0.9 


57±7 


SDP.9 


1 


292.1425 


17.4 


21.7±0.9 


20.3/37.8±2.9 / 


410±50 


11.3 


14.4±1.1 


60±5 


NB.vl.78 


2 


240.2896 


15.4 


36.9±0.4 


5.0±0.8 


510±90 


2.7 


6.7±1.3 


122±24 


G12.v2.30 


2 


231.9458 


4.9 


21.0±0.8 


4.5/4.5±0.7 / 


690±80 


2.7 


6.5±1.2 


128±24 


G15.v2.779 


1 


143.3837 


5.9 


7.3±0.2 




890±140 


3.4 


4.1±0.6 


94±9 


SDP.17b 


2 


298.8646 


28.6 


37.9±1.1 


22.9±2.6 


300±30 


7.2 


7.8±0.9 


85±10 



a: Lines 1 and 2 correspond to the H 2 (2n - 2 02 ) (frest = 752.03 GHz, E up = 137 K) and (2 02 - In) Orest = 987.93 GHz, 
E up = 101 K) transitions, respectively. 

b: zy bs is the frequency corresponding to the zero of the velocity scale of the spectra of Fig. 1. 

c: S^pk corresponds to the peak continuum flux density (mJy/beam); refers to the total spatially integrated continuum flux 
density (mJy). 

d: S^h 2 o pk corresponds to the peak H 2 flux density (mJy/beam); Ih 2 o pk is the peak H 2 velocity-integrated flux density 
(flux in one beam, in Jy kms _1 /beam) and Ih 2 o the spatially integrated H 2 line flux (in Jy kms -1 ). The values derived for 
SDP.81 and G12.v2.30. are underestimated (see text). 

e : /xLh 2 o is the apparent luminosity of the observed H 2 line using the relations given in, e.g., Solomon, Downes & Radford 
(1992). 

f : In the cases of SDP.9, G12.v2.30 and G15.v2.779, the H 2 emission lines are double-peaked and the parameters in the table 
refer to double Gaussian fits. 



due to the low S/N ratio of the H2O data, these numbers are 
less accurate than for most other sources in our sample. Be- 
cause the emission is very extended, the spatially integrated 
H2O line flux Ih 2 o is probably underestimated, as indicated 
by the absence of significant H2O emission seen in the west- 
ern arc, and the very small difference between Ih 2 o and 
In 2 opk compared to the continuum ratio S^/S^pk = 2.2. 

3.2.2. NA.vl.144 at z = 2.202 

The narrow H2O linewidth of NA.vl.144 is the same as is 
seen in the CO emission line (Av ~ 200 kms -1 ). It is the 
narrowest CO linewidth of our sample (Table 3). 

The estimate of the lensing magnification factor from 
detailed modeling based on the SMA image of the 880 /im 
continuum emission indicates a low value of \i ~ 4.6 ± 1.5 
(B13). However, the SMA imaging of this object has lim- 
ited uv coverage (only the extended array was used) which 
might compromise the measurement of the lensing magni- 
fication factor. The corresponding values for the IR and 
CO luminosities appear to be surprisingly high compared 
to the narrow linewidth. Indeed, the uncertain derivation 
of fi from the CO linewidth and luminosity (Harris et al. 
2012) yields /ico ~ 23 ± 15. Therefore, we have doubled 
the uncertainty on the SMA value of fi and adopted the 
rms 2 -weighted average of the estimate from the SMA dust- 
continuum and the CO line, i.e. \i = 5.3±2.9. 

This bright submillimeter galaxy is clearly unresolved 
with our present angular resolution, a result that is con- 
firmed by the small extension observed in the SMA results 
(B13). 

3.2.3. SDP.9 at z = 1.574 

The galaxy SDP.9 is the brightest 250 /im and lowest red- 
shift of the five lenses discussed in Negrello et al. (2010). Its 
lensing magnification factor is estimated to be /i = 8.7±0.7 
from HST imaging (Dye et al. in prep.). We adopted here 



the SMA 880 /im value from B13, /i = 8.5±1.8, which agrees 
with the HST value and also with /ico = 11. 

The H2O intensity of SDP.9 is by far the largest of our 
sample, due to the low redshift of the source. The corre- 
sponding apparent H2O luminosity is among the highest 
found so far. The emission line profiles in H2O and CO 
are similar, including the strong red emission spike and the 
weaker, broader blue emission (Fig. 1). The high angular 
CO data from Cox, Ivison et al. (in prep) as well as the 
HST/WFC3 imaging data (Negrello et al. in prep.) show 
that SDP.9 is a good example of an Einstein ring, with a 
clear asymmetry between the east and west arcs. However, 
our angular resolution is not high enough to resolve the 
emission. 

3.2.4. NB.vl.78 at z = 3.111 

This source has one of the highest redshifts and highestt 
apparent luminosities in our sample. The intermediate 
linewidth of the CO and H2O lines (Av ~ 550 kms -1 ) is 
consistent with those measured in the mid-J CO lines used 
to initially determine the redshift with CARMA (Riechers 
et al. in prep.). The H2O line is indeed the brightest line 
measured in this source so far and thus provides the best 
constraints on the line profile. The moderate CO intensity 
and linewidth indicates intermediate CO, H2O and infrared 
intrinsic luminosities, Lir « 10 13 L . 

The lensing magnification factor, derived from the SMA 
880 /im imaging results, is \i = 10.5±1.4 (B13), a value 
which is adopted in this paper and is consistent with /ico = 
12 (inferred from the CO data of Table 3). 

The source is resolved and shows a relatively complex 
morphology in the H2O and continuum maps. The peak 
value in a single beam, In 2 opk, reported in Table 4, is 
lower than the spatially integrated value Ih 2 o by a fac- 
tor ~2.5. The ratio is similar for the 240 GHz continuum, 
where the peak value is 15.4 mJy/beam and the integrated 
flux density is 36.9 mJy (see Table 4). In addition to the 



southeast main peak, there are northern and southwestern 
extensions (Fig. 2). This morphology is also seen in recent 
SMA 880 /im imaging results (Bussmann et al. in prep.). 
However, the southwest peak is very strong in H2O com- 
pared to the main southeast peak and the peak ratio in 
the continuum. This could suggest that the H2O emitting 
region is more compact than the continuum. 

3.2.5. G12.v2.30 at z = 3.259 

G12.v2.30 is another prominent source among the H- 
ATLAS brightest lenses. Its complex differential lensing has 
been analyzed by Fu et al. (2012), showing that the lensing 
is very different in the near-infrared compared to the sub- 
millimeter or CO emission. The various estimates of the 
lensing magnification factor based on submillimeter images 
or the CO linewidths are comparable: /i = 7.6±1.5 for the 
dust submillimeter and 6.9±1.6forCO(l— 0) imaging by Fu 
et al.; 9.6±0.9 for the same submillimeter data by B13 (the 
value adopted here); /ico = 7±2 by Harris et al. (2012) from 
CO(l — 0) linewidth and luminosity. These values confirm 
that G12.v2.30 is one of the most luminous infrared galaxies 
in H- ATLAS with an estimated L m ~ 1.6 x 10 13 L . 

The H2O profile appears to be double peaked with a to- 
tal width of ~ 700 km/s. Due to the short integration time 
of the observations, the CO results for G12.v2.30 are the 
noisiest of our sample and the details of the profile, in par- 
ticular, its double-peaked nature remain unclear. Within 
the noise, the CO profile is compatible with the H2O emis- 
sion, although additional observations are needed to con- 
firm this result. 

The source is resolved with the present angular resolu- 
tion of 2.0" x 1.0" and shows a complex structure that is 
most clearly seen in the continuum map (Fig. 2): an elon- 
gated arc to the south shows two peaks and emission is seen 
about 2" to the north. Due to the low S/N ratio, only the 
southern component is detected in the H2O emission line. 
G12.v2.30 appears to be a complex source. The dust and 
gas emissions as well as the stellar emission have very differ- 
ent morphologies, which indicates that the lensing potential 
is unusual (Fu et al. 2012). Estimating the total intensity 
of the H2O line remains difficult for G12.v2.30. The peak 
value in a single beam, lH 2 opk, reported in Table 4, is lower 
than the spatially integrated value Ih 2 o by a factor ~2.4. 
The situation is similar for the 240 GHz continuum, where 
the peak value is 4.9 mJy/beam and the integrated flux 
density is 21.0 mJy (see Table 4). This latter value is itself 
much lower than the flux density of ~ 36 mJy measured at 
1.2 mm using MAMBO at the 30-meter telescope (Danner- 
bauer et al. in prep.). However, the reported value for the 
spatially integrated line flux Ih 2 o is less accurate than for 
most other sources in our sample. Because the emission is 
very extended, the spatially integrated H2O line flux Ih 2 o 
is probably underestimated, similarly to the continuum, as 
indicated for the latter by the factor ~1.5 lower compared 
to the MAMBO value. 

3.2.6. G15.v2.779 at z = 4.244 

With a redshift higher than 4 and Fsoo/im = 204±32 mJy, 
G15.v2.779 (also refered to as ID. 141) remains unique in the 
H- ATLAS survey (Cox et al. 2011). Its lensing magnifica- 
tion factor has recently been estimated to be only 4.1 ± 0.2 



(Bussmann et al. 2012), which makes it the most luminous 
H-ATLAS source, with Lir = 2 x 1O 13 L . This is con- 
sistent with the broad CO lines, the widest of our sample 
(760 - 890 kms -1 ), which would imply /i co ~ 8. The H 2 
line has one of the largest S/N thanks to good observing 
conditions in the 2 mm window. 

The width and the characteristics of the complex, 
double-peaked H2O line agree well with those of the CO (7- 
6) and [Cl] lines (Fig. 1). G15.v2.779 has the highest in- 
trinsic H2O luminosity of our sample by a factor ~2. This 
implies that both should have very large H2O luminosities 
if the two components corresponded to two galaxies. 

At the present angular resolution G15.v2.779 remains 
unresolved, in agreement with the compactness of the CO 
(Cox, Ivison et al. in prep.) and the 880 /im continuum 
emissions (Bussmann et al. 2012). 

3.2.7. SDP.17b at z = 2.305 

The detection of the H2O emission line in SDP.17b was 
reported and discussed in Omont et al. (2011). Its H2O line 
profile is in perfect agreement with the CO emission profile 
- see Fig. 1 for the CO(4-3); also Harris et al. (2011) for the 
CO(1-0) and Iono et al. (2012) and Riechers et al. in prep, 
for the CO (3-2) profiles - including the asymmetry with 
the distinct blueshifted emission wing. The source is not 
resolved in the present data and higher angular resolution 
data indicate that the source remains compact and displays 
a velocity gradient (Cox, Ivison et al. in prep.). 

Again the /i value derived by B13, /i = 4.0±0.6, is much 
lower than approximately expected from the CO linewidth 
and luminosity, which yields /ico = 18±8 (Harris et al. 
2012) - note that this latter value is uncertain because of 
the asymmetric CO line profile. However, the image sep- 
arations are at the limit of what can be resolved by the 
SMA. Therefore, we have again doubled the uncertainty on 
the SMA value of /i and adopted the rms 2 -weighted average 
of the estimate from the SMA dust-continuum and the CO 
line, i.e. /i = 4.3±1.2. However, the line of sight toward 
SDP.17b could include an intervening system, as indicated 
by a CO line at z = 0.942 reported by Lupu et al. (2012) 
based on Z-Spec data, in addition to the z = 2.3 line seen 
by the GBT (Harris et al. 2012) and the PdBI (Fig. 1). But 
this CO line at z = 0.942 was not confirmed by Omont et 
al. (2011) and more observations are needed to clarify this 
question. If there was such a dusty object at z = 0.942, the 
lens model derived from the SMA data of this very compact 
source would have to be revised. 

4. Discussion 

4.1. Lensing 

The main purpose of this paper is to investigate the impor- 
tance of the H2O emission in high-z ULIRGs, but not to 
study its spatial structure by separating it from lensing im- 
ages. This is the reason why the observations were carried 
out with the most compact configuration of PdBI to maxi- 
mize the sensitivity for detection. Detailed lensing model- 
ing is thus out of the scope of this paper. However, a good 
control of lensing effects is mandatory for various reasons: 
i) As seen above, the H2O emission could be resolved out in 
sources that are extended, even with the limited resolution 
of the observations reported in this paper, ii) Differential 



Table 5. Properties of the dust continuum and the H2O line emission 



Source ID 


line a obs Lir 


LH 2 0-lm 


LH 2 0-la 


LH 2 0-2m 


LH 2 0-2a 




(1O 12 L ) 


(1O 7 L ) 


(1O 7 L ) 


(1O 7 L ) 


(1O 7 L ) 



SDP.81 


2 


6.1 


(l.sr 


(2.3) 


3.3 


3.3 


NA.vl.144 


1 


11 


9.7 


9.7 


(25) 


(16) 


SDP.9 


1 


5.2 


7.0 


7.0 


(16) 


(10) 


NB.vl.78 


2 


10 


(5.1) 


(7.9) 


12 


12 


G12.v2.30 


2 


16 


(5.8) 


(9.0) 


13 


13 


G15.v2.779 


1 


21 


23 


23 


(55) 


(34) 


SDP.17b 


2 


16 


(8.7) 


(14) 


20 


20 



a: Line 1 is H 2 0(2u - 2 02 ) 752.03GHz and line 2 is H 2 O(2 02 - In) 987.93GHz; 

b: Luminosity values for unobserved lines are written in parenthesis, they are inferred from the other observed lines assuming the 
line intensity ratio from Mrk 231 (Gonzalez- Alfonso et al. 2010, subscript m) or Arp220 (Rangwala et al. 2011, subscript a). As 
detailed in the text, we consider that these values significantly underestimate the H 2 luminosity for SDP.81 and G12.v2.30. 
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Fig. 3. Apparent H 2 luminosity versus apparent infrared 
luminosity of high-redshift sources, for H 2 Line 1 (2n — 2o 2 , 
752.03GHz) with the Arp220 H 2 line ratio (see text), com- 
pared to local ULIRGs (/iLm < 4 xl0 12 L©, Yang et al. in 
prep.). The dashed diagonal lines have constant Lh 2 o /Lir ratios 
6.8 x 10~ 6 and 2.6 x 10~ 6 , equal to those of Arp 220 and Mrk 231. 
Because of their peculiarities, values for SDP.81 and G12.v2.30 
(underestimated by factors possibly up to ^1.5-2 because they 
are very extended, see text), and those of APM 08279+5255 and 
Mrk 231 (type 1 and type 2 QSOs, respectively) are printed in 
blue. 



lensing effects must be controlled for comparison with emis- 
sion from other H2O lines or other molecules or dust. As 
discussed, e.g., by Serjeant (2012), the importance of dif- 
ferential lensing is enhanced for the highest magnifications 
which take place close to caustics, and, especially, for blank 
fields, such as H- ATLAS where the deflectors are either sin- 
gle massive galaxies or small groups of galaxies. A striking 
example is the large difference observed in lensed images 
in the near-infrared compared to the submillimeter or CO 
emission of G12.v2.30 (Fu et al. 2012). iii) Lensing magni- 
fication must be corrected for deriving intrinsic properties 
of the sources such as luminosities or line ratios. 

To infer intrinsic infrared and H2O luminosities of our 
sources, a first approximation is to derive a (mean) lensing 
magnification factor fi. We systematically preferred values 



of fi inferred from submillimeter imaging rather than near- 
infrared HST imaging, since H2O emission should be more 
related to submillimeter dust emission than to near-infrared 
stellar emission, and submillimeter and near-infrared lens- 
ing can be very different as examplified in the case of 
G12.v2.30 (Fu et al. 2012). As discussed in Sec. 3.2, the 
lensing magnification factor /i is thus determined from ex- 
isting lensing models based on 880 jam SMA imaging re- 
sults (B13; see also Bussmann et al. 2012). For two sources 
these values were slightly corrected by also using the values 
of \i derived following a method described in Harris et al. 
(2012), which uses the relation bewteen the CO linewidth 
and luminosity. As the low angular resolution and limited 
S/N ratio of our H2O data do not allow us to trace pos- 
sible differential effects of lensing, we adopted a minimum 
uncertainty of ±20% for all values of fi. 

4.2. Relation between H2O and infrared luminosities 

The intrinsic H2O luminosities, based on the magnification 
factors listed in Table 2, are reported for each source in 
Table 5. The information on the H2O emission is limited 
by the observation of a single emission line per source. This 
prevents us from comparing the H2O excitations between 
different sources. It is nevertheless interesting to compare 
their H2O line luminosity and to check their dependency on 
other properties, such as the infrared intrinsic luminosity, 
Lir, since infrared excitation is thought to play a key role 
in the excitation of the H2O levels. 

In addition, the H2O lines that were observed differ 
from source to source: three sources were measured in the 
H20(2n-2o2) emission line ('Line 1') and four sources in 
the H20(2o2 _ ln) emission line ('Line 2') To deal with this 
disparity, we approximated the integrated line flux, Ih 2 o, 
of the emission line that was not observed, by adopting the 
value of the line ratio of 'Line 2' to 'Line 1' measured either 
in Arp 220, r 2 i =3440/2970 = 1.16 (Rangwala et al. 2011) 
or in Mrk 231, r 2 i = 718/415 = 1.73 (Gonzalez- Alfonso et al. 
2010), listed as 'a' and 'm' in Table 5. The correspond- 
ing values of Lh 2 o are reported in parenthesis in Table 5. 
Finally, the full set of values of the intrinsic luminosities 
Lh 2 o-i/2 f° r lines 1 and 2 and the line ratios of Arp 220 or 
Mrk 231 are reported in Table 5. 
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Fig. 4. Relation of the intrinsic H2O luminosity with the infrared luminosity for the seven ultra- luminous infrared galaxies 
reported in this paper, other well-known high-z sources taken from the literature, and local ULIRGs (Lir < 4 xl0 12 L©, Yang et 
al. in prep.). The luminosity is reported to the H2O (2n — 2o2, 752.03 GHz) emission line (referred to as 'line 1' in this paper), 
where we adopted an Arp220 H2O line ratio when another H2O line ('line 2') was measured (see text). The solid line shows a 
best-fit power-law, Lh 2 o = L?r, where a = 1.17 ± 0.08. Error bars are defined in the following way: i) for Lir, they combine the 
uncertainties on the values of /xLir and the uncertainties on \i given in Table 2; ii) for Lh 2 o, they include the errors on Ih 2 o of 
Table 4, and the same uncertainties on \i except that when the latter is <20% it is replaced by 20% (see Sec. 4.1). Because of 
their peculiarities, values (in blue) for SDP.81 and G12.v2.30 (both underestimated by factors possibly up to ^1.5-2 because they 
are very extended, see text), and those of APM 08279+5255 and Mrk231 (type 1 and type 2 QSOs, respectively) are not taken 
into account for the fit. 



There are several arguments to favor the conversion fac- 
tor of Arp220 rather than Mrk231: for all galaxies, except 
G12.v2.30, the CO excitation seems significantly lower than 
in Mrk 231 (Lupu et al. 2012 ); using the Mrk 231 conversion 
factor seems to yield anomalously low values for the line 1 
luminosity of sources only observed in line 2 (Table 5). 

Therefore, keeping all uncertainties in mind, the rela- 
tionship between Lh 2 o and Lir is displayed in Fig. 4 for 
LH 2 o-ia (line 1 with Arp220 line ratio). This relation- 
ship is inferred from all sources studied in this paper - with 
the exception of SDP.81 and G12.v2.30 whose H2O intensi- 
ties are uncertain and probably underestimated by factors 
possibly up to ~ 1.5-2 because they are very extended - to- 
gether with HLSJ0918 from Combes et al. (2012), and local 
ULIRGs (Lir < 4 xlO 12 L ) from Yang et al. (in prep.). A 
best-fit power-law to this relation following 

£h 2 o = Lf R (1) 

yields a = 1.17±0.0.10 using 'line 1' and the Arp 220 
line ratio. Using the Mrk 231 line ratio would 
yield a = l.ll±0.10, whereas using 'line 2' would give 



a = 1.22±0.10 and 1.28±0.12 for the Arp 220 and Mrk 231 
line ratios, respectively. Independent of the selected lines 
or ratios, these results indicate that Lh 2 o increases rapidly 
with Lir, possibly faster than linear. We note, however, 
that since systematic errors were not taken into account in 
this analysis, adding the statistical and systematic errors 
could result in a relation that is closer to being linear. A 
steep increase is consistent with a significant role of infrared 
radiation in the excitation of the water vapor lines. It is in- 
teresting to note that a similar behavior with the infrared 
luminosity has been found for the intensity of local H2O 
and OH megamasers, which could be interpreted as a re- 
sult of the importance of infrared pumping on these masers 
(e.g. Lo 2005). 

However, it is important to stress that, as shown e.g. 
by the modeling of H2O excitation in Mrk 231 by Gonzalez- 
Alfonso et al., the intensity of (optically thick) H2O lines 
results from a complex interplay between various parame- 
ters that includes gas temperature and density, H2O abun- 
dance, infrared radiation field, spatial distribution of H2O 
and dust, etc. The present data, with only one H2O emis- 
sion line per source, prevents us for constraining the water 



excitation. Multi-line H2O observations are required for a 
detailed modeling of the excitation. Therefore, we defer the 
study of the water excitation in these high-z sources to a 
future paper, which will report measurements of multiple 
water transitions. 

The high-J HCN lines seem to display a similar behav- 
ior with the infrared luminosity, with Lhcn growing possi- 
bly faster than linear with Lir. For instance, Bussmann 
et al. (2008) reported that L H cn(3-2) varies as Ljj^ ±01 
in a sample of 30 nearby galaxies, although correcting 
for the difference in emission extension seems compatible 
with a linear relation (Gao et al. in prep.). The strong 
HCN(5 — 4) emission observed by Weiss et al. (2005) in 
APM 08279+5255 seems to require infrared excitation. The 
high HCN 6-5/5-4 ratio observed by Riechers et al. (2010) 
ultimately confirmed that infrared excitation is important 
for the HCN excitation in this source. Note that this is 
not true for HCN(1 — 0), which gives a strictly linear rela- 
tion, with an index of 1.0, consistent with the interpretation 
that HCN(1 — 0) is proportional to the mass of dense gas 
(Solomon et al. 1992; Gao & Solomon 2004). However, 
because the H2O and high-J HCN line intensities may be 
increasing faster than linear with Lir may indicate the im- 
portance of IR excitation. This needs to be confirmed by 
better statistics and by future observations of higher tran- 
sitions. 

4.3. Comparison of H 2 lines with high-J CO lines and the 
continuum 

In all cases where the information is available on adjacent 
CO lines - CO(6-5) or (7-6) for H 2 O(2ii-2 02 ) (line 1), and 
CO(8-7) or (9-8) for H 2 O(2 02 -lii) (line 2) - the line flux 
IH20-1/2 is a ^ least ~0.3-0.5 that of the neighboring high- 
J CO lines. The lack of detailed data on high-J CO emis- 
sion for most of the sources limits any detailed comparative 
study. 

The three SDP sources and G12.v2.30 benefit from a 
complete spectral coverage from 200 to 300 GHz with ZSpec 
(Lupu et al. 2012 and in prep.). This range includes 2 
to 5 high-J CO lines, depending on the redshift. Even 
though the sensitivity and the spectral resolution of ZSpec 
are limited, these data allow us to derive useful trends on 
the strength of CO as the approximate spectral line energy 
distributions (SLEDs) derived by Lupu et al. (2012). 

In SDP.17b, the ratio I(H 2 O(2 02 - ln))/I(CO(8 - 7)) 
is about 0.5 (Omont et al. 2011). For SDP.9, informa- 
tion is lacking on CO emission from levels J > 6 be- 
cause of the low redshift (z = 1.57) of the source. How- 
ever, the CO (6-5) emission line is strong with a ratio 
I(H 2 0(2n -2o 2 ))/I(CO(6-5))~0.36. In the case of SDP.81, 
the H 2 and CO emission lines are both weak. Accounting 
for the part of the H 2 emission that is not detected in the 
intrinsically weaker western arc (see Sect. 3.2.1), we esti- 
mate that the ratio I(H 2 O(2 02 - ln))/I(CO(7-6)) is w 0.3 
- note that the CO J > 7 lines are not detected in the ZSpec 
data reported by Lupu et al. (2012). The ZSpec spectrum 
of G12.v2.30 displays strong lines of H 2 0(3 2 i — 3i 2 ) and 
CO (10-9), with about similar intensities (Lupu et al. in 
prep.), but their ratio may be uncertain by a factor up to 
two because of the noise in this frequency range of ZSpec. 

The best-documented comparison of H 2 and CO 
lines is found for G15.v2.779 because of its high red- 
shift (z = 4.24) and the sensitive PdBI CO observa- 



tions of Cox et al. (2011). For G15.v2.779, the ratio 
I(H 2 0(2n - 2 02 ))/I(CO(7 - 6) = 0.52. For the two other 
sources, NA.vl.144 and NB.vl.78, the lack of information 
on the high-J CO levels prevents any comparison with the 
H 2 emission lines. 

Moreover, the intensities of H 2 lines (Ih20 - Table 4) 
and of the lower- J CO 2 mm and 3 mm lines (Table 3 and 
Fig. 1) are roughly comparable. Another direct comparison 
is the relation between the intensity of the H 2 emission 
line (hence the H 2 luminosity) with the underlying con- 
tinuum emission (hence the infrared luminosity) - see Table 
4 and Figs. 3 and 4. But the comparison between Ih 2 o and 
the continuum flux density S u significantly depends on the 
redshift and the rest frequency. Finally, as shown by Fig. 
1 and Table 4, the H 2 and continuum flux densities, e.g. 
Sz,H 2 opk and S^pk, are roughly comparable, although their 
ratio obviously depends on the linewidth. 

4.4. Implication of ubiquitous detection of H 2 in high-z 
ultra-luminous infrared galaxies 

As in the case of SDP.17b (Omont et al. 2011), the strength 
of the H 2 2o 2 -ln and 2n-2o2 emission lines in six new 
high-z ultra-luminous galaxies implies that water vapor has 
a high abundance in these sources - probably n(H 2 0)/n(H 2 ) 
> 10 -6 as in Mrk231 (Gonzalez- Alfonso et al. 2010), al- 
though this has to be confirmed by detailed modeling - and 
that its rotation level ladder is excited up to energies of at 
least ~ 100 K. This is consistent with the recent results 
obtained on local ULIRGs such as Arp220 and Mrk231 
(Rangwala et al. 2011; van der Werf et al. 2010; van der 
Werf et al. in prep.; Yang et al. in prep.). 

Our results seem to indicate that high-z ultra-luminous 
galaxies are analogs of their local counterparts scaled by 
factors up to an order of magnitude in infrared luminos- 
ity and star formation rate. The similitude in excitation 
of high-energy H 2 levels, observed in local ULIRGs or in 
APM 08279+5255 (van der Werf et al. 2011), still needs 
to be confirmed in high-z sources by future observations 
of the higher frequency H 2 lines connected to these lev- 
els. Nevertheless, the results reported here already indicate 
many similarities, such as the possibly slightly faster-than- 
linear increase of the intrinsic H 2 luminosity, Lh 2 o, with 
the infrared luminosity Lir; the corresponding fact that 
the ratio Lh 2 o/Lir in high-z sources is slightly higher than 
in local ULIRGs (Fig. 3); and the high intensity ratio of 
H 2 lines to nearby high-J CO lines, in between 0.3 to 
0.5, as observed in local ULIRGs. It is therefore likely that 
the infrared radiation plays a significant role in the H 2 
excitation, as demonstrated in the case of Mrk231 by the 
detailed modeling of its spectrum by Gonzalez- Alfonso et 
al. (2010). 

The ubiquity of H 2 emission lines in the submillimeter 
spectra of luminous high-z starburst galaxies underlines the 
importance of water vapor in probing their warm, dense, 
and dusty interstellar cores and, in addition, their strong 
local infrared radiation fields. Indeed, the H 2 lines con- 
stitute a totally different diagnostic from, e.g., CO or low-J 
HCN emission lines, that are excited by collisions and hence 
trace only temperature and density. Because of the large 
electric dipole of H 2 and its large rotational constant, the 
critical density for collisional excitation of its excited lev- 
els is extremely high, typically 10 8 cm -3 , much higher than 
expected densities in our sources. The relatively easy de- 



tection of H2O emission lines in the sample of high-z star- 
burst galaxies studied in this paper is therefore indicative 
of the presence of dusty, infrared-opaque nuclear disks in 
their centers as in local ULIRGs and of a significant role 
of infrared excitation of H2O (e.g. Gonzalez- Alfonso et al. 
2010). 

The role of a strong AGN in the excitation of H2O of 
local ULIRGs is not yet entirely clear (van der Werf et al. 
2010). There is no definite evidence of the presence of a 
strong AGN in any of the seven sources discussed in this 
paper, neither from their optical/near- infrared spectra nor 
from their X-ray emission. There is also no indication of 
a high radio or mid-infrared excess that could have been 
detected in the FIRST radio survey (Becker et al. 1995) or 
the Wide-field Infrared Survey Explorer (WISE; Wright et 
al. 2010) - see Omont et al. (2011) for a detailed discussion 
of the radio emission of SDP.17b. However, in the case of 
SDP.81, the radio spectral index suggests that a fraction of 
the radio emission is powered by an AGN (Valtchanov et 
al. 2011). 

We stress again that most of our sources are much more 
powerful than local ULIRGs and that they probe a regime 
of infrared luminosities up to a few 10 13 L , i.e., an order 
of magnitude higher than local ULIRGs. It is probable that 
the local radiation field is also stronger, and that it might 
reach the limit of stability of such cores against radiation 
pressure, possibly close to the limit of maximum starbursts. 
As starbursts with intrinsic far-infrared luminosities much 
higher than a few 10 13 L do not appear to exist (see e.g. 
Karim et al. 2012), it seems that some of the objects that 
we are studying, such as G15.v2.779, might be close to this 
critical regime. 

All sources studied in this paper are highly magnified 
by gravitational lensing. This makes the detection of lines 
such as H2O much easier, especially for the highest lens- 
ing magnification factors, i.e., the lowest intrinsic luminosi- 
ties. However, for the strongest luminosities, a few 10 13 
L , lensing magnification factors go down to a few units; 
e.g. fi ~ 4 for G15.v2.779. Therefore, for similarly lumi- 
nous sources, H2O lines should be easily detectable even in 
the absence of lensing magnification. For instance, a source 
identical to G15.v2.779 (L m - 2 x 10 13 L at z > 4) would 
have Ih20 ~ 1 Jy.km/s without any magnification, which is 
well within the reach of the PdBI in its current configura- 
tion. 



5. Conclusion 

This paper reports a significant step forward in the study 
of water vapor emission in the warm and dense cores of 
ultra-luminous starburst galaxies at high redshift. The 
seven sources studied here are strongly lensed galaxies dis- 
covered in the Herschel H-ATLAS survey and, using the 
PdBI, all were detected in H2O either in the (2o2 _ ln) or 
the (2n-2o2) rotational line. One of the sources (SDP.17b) 
was previously reported in Omont et al. (2011). In all 
cases, the H2O emission lines are strong, with integrated 
line fluxes ranging from 1.8 to 14Jykms _1 , comparable 
both in strength and in the characteristics of their profiles 
to the CO emission lines. The derived apparent infrared 
luminosities, fiL^o, range from ~3 to 12 x 1O 8 L . WhenD 
correcting for the magnification factors /1, which are esti- W 
mated from lensing models based on submillimeter imag- [5] 
ing and corroborated by the empirical relationship between 



CO luminosity and linewidth, the derived H2O luminosi- 
ties show a strong dependence on the infrared luminosity. 
It is found that Lh 2 o varies as ~Lj^?, the exact shape of 
the relation remaining uncertain mainly due to the uncer- 
tainties in the derivation of ja. This slightly non-linear re- 
lation may indicate that infrared pumping plays a role in 
the excitation of H2O lines. Water emission is therefore 
expected to become very strong in the most ultra-luminous 
(Lir > 10 13 L ) galaxies at high redshift. 

The results described in this paper underline the fact 
that the most powerful high-z infrared galaxies are analogs 
of local ULIRGs scaled by factors up to 10 in infrared lumi- 
nosity and star formation rate. The fact that the excitation 
of H2O levels with E up ~ 100-150 K in high-redshift lumi- 
nous galaxies appears to be similar to what is seen in local 
ULIRGs is supported by the present results, including the 
rapid increase of Lh 2 o with Lir, the slightly higher ratio 
Lh 2 o/Lir in high-z sources than in local ULIRGs, or the 
high intensity ratio of H2O lines to nearby CO lines. How- 
ever, future observations of the higher frequency H2O lines 
connected to upper levels in high- z galaxies are needed to 
compare the excitation of the high-energy H2O levels with 
E up > 150-200 K to that of local ULIRGs. 

The detection of H2O emission lines in seven high- 
redshift luminous galaxies reported in this paper demon- 
strates that the millimeter /submillimeter lines of water 
vapor are a new and key diagnostic tool to probe both 
the extreme physical conditions and the infrared radiation 
in the warm dense environments of luminous high-z star- 
burst galaxies. Follow-up observations of such sources, tak- 
ing advantage of the increased sensitivities that are and 
will become available with ALMA and the upgraded PdBI 
(NOEMA) facilities, will allow more comprehensive studies 
of water vapor in high-z luminous galaxies, including un- 
lensed candidates, providing essential and new insights into 
the physical conditions of these galaxies which were present 
when the universe was young and have no equivalent in the 
local universe. 
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